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Abstract 

A global analysis of the solar neutrino data from all solar neutrino ex- 
periments combined with the KamLAND data is presented assuming that 
the solar neutrino deficit is due to the matter-enhanced spin-flavor preces- 
sion effect. We used two types of magnetic field profiles throughout the 
entire Sun: Wood-Saxon shape and the Gaussian shape. We showed that 
for Dirac neutrinos, the allowed regions are independent of the magnetic 
field profiles for all of the magnetic moments that we used in this paper 
and the allowed region in the large mixing angle (LMA) region shifted to 
the small mixing angle (SMA) region as fiB value is increased. We calcu- 
lated the allowed regions at 95% CL. We also find a limit for the electron 
neutrino magnetic moment at 0.95CL so that {iB < 0.2 x 10" ' fJ, B G for 
both magnetic field profiles at lcdevel. 

1. Introduction 

The earlier solar neutrino experiments (Homestake and all three Gallium 
experiments) [1 — 4] showed deficits of the neutrino flux from the Sun when com- 
pared to the standard solar model(SSM) predictions [5]. 

One of the most popular solutions of this problem is MSW effect [6, 7]. In 
this solution, when electron neutrinos pass through the electronic matter, they 
are exposed to matter effects which can cause an almost complete conversion of 
electron neutrinos to another neutrino types. 

Another solution is the magnetic moment solution which we investigated in 
this paper. If the neutrinos have large magnetic moments and pass through a 
region with a magnetic field, the helicity of the neutrino can be flipped, then this 
spin flip changes the left-handed electron neutrino to a right handed electron 
neutrino that we cannot detect [8 — 11]. 

Lim and Marciano examined [12] the combined effect of matter and magnetic 
fields on neutrino spin and flavor precession. They developed the idea of the 



resonance spin flavor precession(RSFP) in 1988. In this solution based on the 
RSFP mechanism, the combined effect changes the neutrino's chirality and fla- 
vor simultaneously. Matter-enhanced spin-flavor precession of solar neutrinos 
with transition magnetic moments for chlorine and gallium experiments was in- 
vestigated in detail by Baha et al [13]. In recent years there have been several 
other studies on RSFP investigating different aspects[14 — 23]. Chauhan et.al 
[23] examined the combined action of neutrino oscillation and spin-flavor preces- 
sion (SFP). They considered the flux of electron antincutrinos coming from the 
Sun as a possible observable effect of SFP, and as a result tried to put an upper 
bound as \xB < 2.8 x 10~ 19 MeV. Also to examine in detail the effect of the 
Sun's magnetic field, a statistical analysis of the solar data has been performed 
in [21] and from the minimum of chi-squares, the value of the magnetic field for 
different profiles can be deduced. In addition to these, one can examine first 
the allowed region of neutrino parameter space for each neutrino experiment 
(Chlorine, Gallium, SK and SNO) depending on the magnetic field strength in 
the Sun and secondly can find the allowed regions from combined solar v exper- 
iments at the same magnetic field strength, by adding the KamLAND data to 
the combined solar v data to make the analysis global. In our present present 
work we followed such a path to extract a upper /iB value. 

In the present work a global analysis of the solar neutrino data combined 
with the new binned KamLAND data [24] is presented assuming the solar neu- 
trino deficit to be resolved by the matter-enhanced spin-flavor precession. We 
calculated survival probability for two different magnetic field profiles. 

Standard least-squares analysis of solar neutrino data is often used to obtain 
information on the values of the allowed regions for the oscillation parameters 
namely Am 2 , tan 2 9. Large mixing angle (LMA), small mixing angle (SMA), 
and low Sm 2 (LOW) regions are usually known as MSW solutions. Besides 
the earlier solar neutrino experiments Homestake [1] and Gallium[2 — 4], SNO 
charged-current and neutral current [25, 26] and Super Kamiokande [27] results 
comfirmed solar neutrino deficit. 

In 2002 global analysis of all solar neutrino experiments showed that the 
LMA solution was the most likely solution in the neutrino parameter space [28]. 
Also data from the KamLAND reactor neutrino experiment indicated the LMA 
region [24, 29]. So, solar and reactor neutrino experiments strongly pointed out 
the LMA region. 

In our paper, we examined how the solar magnetic field in the Sun can 
change this region. We found results on the allowed regions and chi-square 
values for two different magnetic field profiles and for nine different values of 
neutrino magnetic moments. Our results showed that this allowed region in the 
LMA shifted to the SMA region for certain ^.B values. From these observations 
we were able to put upper limits on the ^iB values. In section 2, we give general 
information about equations that governs the neutrino propagation assuming 
that neutrinos have magnetic moments. Properties of the magnetic field profiles 
are given in section 3. We give detailed statistical analysis in section 4. Finally 
our results and conclusion are given in the section 5. 

2. Matter-enhanced spin-flavor precession in the Sun 
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For the electron neutrino, in the case of v PL — ► v eR precession the v eL and 
v ea interact differently with matter. The differences in their interactions with 
matter suppresses precession by splitting their degeneracy. 

The evolution equation that describes the propagation throughout matter of 
the two chiral components v eL and v PR with magnetic moment fj, Bv is [8,9,30] 



dt 







' v e (t) 


/jlB ' 




















(1) 



where B is the transverse magnetic field and V e (t) is the "matter" potential 
that is the contribution of matter to the effective mass. In the standard model 
for an unpolarized neutral medium, 



V e (t) = ^(2N e -N n 



(2) 



where N e and N n are electron and neutron number densities respectively and the 
Gf — 1.16636 x lO~ 5 GeV~ 2 . Also in the Sun electron and neutron number den- 
sities are well approximated by N e ~ 6N n ~ 2.4 x 10 26 exp(— r /Q.Q%Rq) / cm? . 
In the [30], the resonance condition was found as N e ~ N n /2. This neutron 
number density condition can not be found within the Sun, but it can exist in 
a supernonova [31,32]. 

Although there is no resonance region in the Sun for spin-precession, it was 
shown that there is a resonance region where neutrino spin-flavor precession 
may occur for a medium with changing density. 

To explain the resonant spin-flavor precession, we first consider two genera- 
tions. For definiteness, we examine v e — system. The evolution equation for 
a neutrino that propagates through matter and a transverse magnetic field B is 



dt 
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where 2x2 submatrices are 
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where 9 is the mixing angle, Am 2 is the difference of the mass and E is the 
neutrino energy. The matter potentials for a neutral unpolarized medium are 
given as 



Ve(t) 



V.. 



2l 

V2 



(6) 
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Figure 1: Magnetic field profiles: (a) Wood-Saxon shape; (b) Gaussian shape 



3. Magnetic Field Profiles 

In our analysis, we used two types of magnetic profiles. 
First, we took the magnetic field profile to be a Wood-Saxon shape of the 
form 



B(r) 



Bn 



l + exp[lO(r-i? )/i? ] 



(7) 



where B is the strength of the magnetic field at the center of the Sun. 

The second magnetic field profile that we used is a Gaussian shape. Altough 
there are many other profiles [14, 17 — 19, 21, 23] we considered only those two as 
typical for our purpose as shown figure 1. The effects of these different profiles 
have been discussed in the conclusion section. 

4. Statistical Analysis 

In the literature, there is a common way often called \ 2 analysis to find the 
values of the neutrino oscillation parameters Am 2 , tan 2 9 and to calculate the 
confidence levels of allowed regions and the goodness of a fit [33 — 36] . In our 
analysis, we use "covariance approach" to find the allowed regions mentioned 
above. By this method, one minimizes the least-squares function 



(oxp) _ jy(thr)^ 



(8) 



where 1 is the inverse of the covariance matrix of experimental and theoretical 
uncertainties, i^ exp ) is event rate calculated in the ith experiment and 
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is the theoretical event rate for ith experiment. The indices indicate the solar 
neutrino experiments: = 1, ...,N exp with N cxp = 4. 

V ili2 = V ili2 (exp) + V ni2 (thr) (9) 

V ili2 (exp) = ^ li2 < xp < xp (10) 
where c° xp are the experimental uncertainties for ith experiment. 



V ili2 (thr)=5 ili2 ^R^ )2 ^C%f + (11) 
ii=i 

E <^ )2 < 2 t )2 E^ fe ^ fe (AlnX fc ) 2 
iu'2=i fc=i 

where the indices jij'2 = 1, ...,8 indicate solar neutrino fluxes produced in the 
eight thermonuclear reactions in the Sun: pp, pep, hep, 7 Be, 8 B, 13 N, 15 0, 17 F, 
respectively. The index k=l,...,12 denotes the input astrophysical parameters 
X k in the Standard Solar Model(SSM), on which SSM neutrino fluxes $| SM 
depend. The logarithmic derivatives 



91n$f SM 

govern the uncertainties of the neutrino fluxes <j>? 5 ' M . AlnX^ and AlnCj*' 1 '*'' 
are la relative uncertainties of SSM input parameters and the energy-averaged 
cross section( C^^), respectively. 

Theoretical event rates for the radio-chemical experiments, chlorine exper- 
iments(Homestake) and gallium cxperimcnts(SAGE,GALLEX,GNO), can be 
found as 

Rf r) = J dEfr (E) aj^Pij ( Ve (13) 



such that 



R (thr) = R fhr) (w) 
i=l 



where fa (E) is the flux at energy E coming from ith reaction and <jj (E) is the 
cross section for detector j . 

Solar neutrinos are observed in SK via the Cerenkov light from the neutrino- 
electron scattering (ES) reaction: 

v x + e~ — > u x + e~ 
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where v x can be v e or v^. 

Since SNO is heavy water-Cerenkov detector, it observes solar neutrinos by 
charged-current(CC) and neutral-current (NC) in addition to ES 



v e +d^p + p + e- (CC) 
v x +d^n+p + v x (NC) 

The Cerenkov light is generated by recoiling electron from the ES and CC 
reactions. 

Because of higher threshold energy of SK and SNO experiments(T t / (res ^ jrf > 
5 MeV), they are sensitive to only 8 B and hep neutrinos. Due to the small flux 
of hep neutrinos we completely neglect its contribution to the total rates. 

For SK and SNO, the theoretical event rates from ES: 

R E S S K = J dEfr (E) {ai e (E)P j (v e ^u e ,E, t)+ (15) 
ai x (E)[l-Pj(iy e ^u e ,E,t)]} 

where 

<^ {E) 'L nm iT < 16 > 

Here T is the kinetic energy of the recoiling electron. T m ; n = E t hreshoid — m e 
with Ethreshoid = 5.5 MeV and Tmax — IE / (IE + Trie) are the minimum and 
maximum kinetic energy of the recoiling electron, respectively. The differantial 
cross sections for v e — e~ and v x — e~ scatterings are [37] 

<]1 " °elgi + .4(1 - T/Ef g L g R (T/E 2 )] (17) 



dTdE 
with 

9L = (±^ + sin 2 w ) gR^sii^ew 

The upper sign is for v e — e~ scattering and the lower sign for v x — e~ scattering; 
x can be /i or t. The cross section factor is 

(7 e - = 88.083 x 10- 46 cm 2 

For SNO, in addition to the ES, the theoretical event rates come from the 
CC and NC reactions 

R CC = J dE<j>i (E)^ cc (E)P J {y e - v e ,E,t) (18) 
R NC ' = J dE^(E)^ NC {E) (19) 
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Total event rates are 



ftSNO _ ftCC _|_ ftNC 

Because of its substantially lower cross section compered to the other two reac- 
tions', we don't take into account R ES . 

We took fluxes and cross sections for event rates and error matrixes from 
Bahcall [37, homepage]. 

We need to calculate XKamLAND f° r tne global analysis 

(20) 

For this purpose we note that, the goal of the KamLAND is to search for 
antineutrinos emitted from distant power reactor through the reaction 

p + T7 e ->n + e+ (21) 
As in [38] , the antineutrinos' energy spectrum is given by 

dNi 

— ^ oc exp(a a + ai E + a 2 E 2 ) (22) 
dhj 

where j = 1,2,3,4 corresponding to the four isotopes 235 U, 239 P, 238 U, 241 Pu 
and the fitted values of the a& are given in detail in [39] for each isotope. 

In the case of two generations, the survival probability for the electron an- 
tineutrinos coming from the ] th reactor 

. 2™ . 2 /l-27Am 2 ( e y 2 )d ? (A:m), 
P(77 e ^77 e ) = l-sin 2 20sin 2 ( ^QeV) } ( } 

here dj is the reactor-detector distance. The number of expected events for each 
bin at KamLAND is 



Nf hr (t,E, Am 2 , sin 2 29) = V N p J dE v j dE e R(E v ,E e )x (24) 



E 



a{E v )P{V e ^V e ,E v ) 



. 47Td 2 

J 3 

here N p is the number of free protons in the fiducial volume of detector and 
r\ is the efficiency which are given in [24]. Sj is the initial energy spectrum 
of reactor j calculated using the thermal power and the isotropic composition 
of each detector that is given in detailed in [39] and (r(i?^) is the lowest cross 
section 

a{E v ) = -^^ Pe E e (25) 

where / = 1.69 is the integrated Fermi function for neutron, m e is the positron 
mass, r„ is the neutron lifetime, p e and E e are the positron momentum and 
energy, respectively. So that we have 
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E e = E v - 1.293 MeV 



(26) 



The energy resolution function R(Ev,E e ) which depends on visible energy 
(Ey) and true positron energy (E e ) is given by 

1 (E v -E c + mc ) 2 

R(E v ,E e ) = ^==e (27) 
V 27rcr o 

in which E v = E e + m e and er = 6.2%a/BJ [24]. 

The new KamLAND data results were reported at [24] in 13 bins above the 
threshold which is 2.6 MeV. To find x 2 for KamLAND spectrum data, due to 
the low statistics, we use x 2 assuming a Poisson distribution given by 



/V exp ' 



+ ^ (28) 



where sum is over the KamLAND spectral bins, G 1 sys is the systematic uncer- 
tainty taken to be 6.5% and k is allowed to vary freely. 
5. Results and conclusions 

In our calculations, we assumed that the magnetic field extends over the 
entire Sun for either Wood-Saxon shape or Gaussian shape. We use the neutrino 
spectra from the Standard Solar Model of Bahcall and his collaborators [40] . In 
our analysis we calculated allowed regions at 95% confidence level. 

We first considered only solar neutrino data. In our statistical analysis, 
we used the covariance approach. We show the allowed regions of neutrino 
parameter space in figure 2 for each solar neutrino experiment seperately at 
four /iB values using Wood-Saxon field profiles. In figure 2 each column and 
row are for the same experiment and at the same [xB value respectively(e.g. in 
the second row at the third column, an allowed region for SK experiment at 
fiB = 2 x 10~ 7 yU B G is seen). 

In figure 3, we displayed the allowed regions from combined solar neutrino 
experiments at the same values of \iB in figure 2. 

After we investigated the combined solar data at different \iB values, we 
looked for allowed regions at different confidence levels from new binned Kam- 
LAND data [24] in figure 4. We show the allowed regions from our global 
analysis combining solar and new KamLAND data in figure 5 for nine differ- 
ent \iB values(0, 1, 1.5, 2, 2.5, 3, 4, 5, 10 x 10~ 7 fi B G ). In that figure, our results 
showed that the allowed regions in the LMA shifted to the SMA region as iiB 
value is increased. This shift enables us to put an upper limit on the value of 
\iB , since the latest experimental data prominently indicate the LMA region 
(namely excluding SMA and the others). 

Tablel. Best fit points of global analysis for both magnetic field profiles. 
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1.77xl0~ 4 
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42.45 
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1.77xl0~ 4 


0.057 


49.39 
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1.53xl0~ 4 


0.022 


53.78 


54.13 


10 


1.53xl0~ 4 


0.012 


55.5 
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We next calculated the global analysis again for a magnetic profile of Gaussian- 
shape to show how the magnetic profile in the sun effects the allowed regions 
and the minimum chi-squares. The results of our global analysis for both mag- 
netic field profiles are given in table 1. As can be seen from the table there are 
no appreciable differences between the effects of the two magnetic field profiles 
on the allowed regions and chi-squares values. 

Finally in figure 6 we present projection of the global A\ 2 function on the 
fiB to place a limit on /j,B for the two magnetic field profiles. From this 
figure we found a general limit, since the graphes for the two magnetic profiles 
almost coincide: [iB < 0.2 x 10~ 7 [i B G, 0.5 x l0~ 7 [i B G, 1.0 x 10~ 7 [i B G for 
la, 2a, 3tr limit, respectively. On the other hand there is intensive work related 
to the electromagnetic properties of the neutrinos. Since neutrinos are known 
to be massive they can get tiny magnetic moments even in the SM also in 
its extensions. Astrophysical considerations [41] put strong constraints on the 
magnetic moment of the neutrinos as \x v < 10~ 12 fi B , however these bounds are 
model dependent. Also the reactor experiments bring less restrictive limits of 
/i„ < 1.3xlO~ 10 /i B [42,43]. An analysis of the Super-Kamiokande solar neutrino 
data [44] enables the authors to put again a similar limit of \i v < 1.5 x 10~ 10 /i B . 
As for the magnetic field of the Sun, from the existing observations an upper 
limit of 10 7 G at the core and a maximum magnitude of 3 x 10 5 G at the bottom 
of the convective zone are usually taken in the literature [21, 22]. Hence the 
above restriction according to our analysis on the /iB value agrees with the 
present limits of \i v and Bg un . 

Also the consideration of the two magnetic field profiles seem to be satis- 
factory at this level, because there are no differences between their effects on 
the allowed regions and the chi-square values for both of them. Actually the 
whole analysis can be repeated for many other field profiles [14, 17 — 19, 21, 23] 
in order to better investigate the detailed behaviour of the magnetic field in the 
Sun. This could be a topic of a future work. 
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tan 2 (6) tan 2 (9) tan 2 (9) tan 2 (9) 

Figure 2: The allowed regions of neutrino parameter space for each solar 
neutrino experiment seperately at fiB = 0, 2, 5, 10 x 10~ 7 [i B G and at 95% 
CL. Each column and row arc for the same experiment and at the same \iB 
value, respectively(e.g. in the second row at third column, an allowed region for 
SK experiment at fiB = 2x 10~ 7 /x B G is seen) 
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Figure 3: The allowed regions from combined solar neutrino experiments ( chlo- 
rine, all three gallium, SK and SNO experiments) at the same [iB values and 
CL as in figure 2. Stars indicate the local best-fit points. 
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Figure 4: Allowed regions from the new KamLAND spectrum at different con- 
fidence levels. The star indicates the best fit point. 
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Figure 5: Allowed regions from the combined solar+new KamLAND spectral 
analysis at different \iB values and at 95% CL. Stars indicate the best fit points. 
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Figure 6: Projection of the global A% 2 function on /j,B. The solid line and short- 
dashed lines represent the calculation using Wood-Saxon shape of magnetic 
profile and Gaussian shape of magnetic profile, respectively. 
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